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Experiments on high-energy heavy-ion collisions reveal the formation and some intriguing properties of the 
Quark-Gluon Plasma (QGP), a new phase of matter predicted by Quantum Chromodynamics (QCD), the quan- 
tum field theory of strong interactions. The phenomenological success of relativistic hydrodynamic simulations 
with remarkably weak shear viscosity, modeling QGP as an almost perfect fluid, are in favor of the occurrence 
of a strongly-coupled QGP expanding and cooling during the reaction. A derivation of these features in QCD at 
strong coupling is still lacking and represents a very intricate theoretical challenge. As a quite unique modern tool 
to relate these dynamical features to a microscopic gauge field theory at strong coupling, time-dependent real- 
izations of the AdS/CFT correspondence provide a fruitful way to study these properties in a realistic kinematic 
configuration. Relating a 4-dimensional Yang-Mills gauge theory with four supersymmetries (which is a conformal 
field theory, CFT4) with gravity in Anti-de Sitter space in five dimensions (AdSs), the AdS/CFT correspondence 
provides a useful "laboratory" to study yet unknown strong coupling properties of QCD. Besides the interest of 
revealing new aspects of the AdS/CFT correspondence in a dynamical set-up, the application to plasma formation 
leads to non trivial theoretical properties, as we will discuss in the lectures. The highlights of the present lectures 
are: 

1. Emergence of an (almost) perfect hydrodynamic fluid at late proper-times after the collision. 

2. Duality between an expanding 4-dimensional plasma and a black hole moving radially in the bulk. 

3. Intimate link between conformal hydrodynamics and Einstein's equations in the asymptotically AdSs space. 

4. Possibility of studying the far-from-equilibrium stage of a gauge field theory at early collisional proper-times. 



1. Heavy-Ion reactions and the QGP 

One of the most striking lessons one may draw 
[112] from experiments on heavy-ion collisions at 
high energy (e.g. at the RHIC accelerator) is that 
fluid hydrodynamics seems to be relevant for un- 
derstanding the dynamics of the reaction (see, for 
instance, the reviews |5]). Indeed, the elliptic flow 
[3] describing the anisotropy of the low transverse 
momentum particles produced in a collision at 
non zero impact parameter implies the existence 
of a collective flow of particles. It agrees with 
the picture of an hydrodynamical pressure gra- 
dient due to the initial eccentricity in the colli- 
sio Moreover, the hydrodynamic simulations 
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3 Note that a preliminary analysis of recent LHC data on 
the elliptic flow at the Alice detector seems to confirm the 



which are successful to describe this elliptic flow 
are consistent with an almost "perfect fluid" be- 
havior, i.e. a small viscosity over entropy ratio r\j s 
[5]. This ratio is particularly interesting since it 
depends mainly on fundamental features of the 
fluid, being in particular independent of the par- 
ticle density (see later, eq.([TJ). 

The validity of an hydrodynamical description 
assuming a quasi-perfect fluid behavior has been 
nicely anticipated by Bjorker@ in Ref.[5|. The so- 
called Bjorken flow is based on the approxima- 



validity of the hydrodynamic approach [2]. 
4 The introduction of hydrodynamics in the description of 
high-energy hadronic collisions has been originally pro- 
posed by Landau [6] , assuming "full stopping" initial con- 
ditions which result in a non boost-invariant solution or 
Landau flow (see [7] for a unified description of Bjorken 
and Landau flows). However "full stopping" initial condi- 
tions do not seem to agree with present day data. 
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Figure 1. Description of QGP formation in heavy 
ion collisions. The kinematic landscape is defined 
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The coordinates along the light-cone are xo ±a;i, the 
transverse ones are {x2, £3} and r is the proper time, 
while rj is called the "space-time rapidity" . 



tion of an intermediate stage of the reaction pro- 
cess consisting in a boost-invariant quark-gluon 
plasma (QGP) behaving as a relativistically ex- 
panding fluid. It is formed after a (strikingly 
rapid) thermalization period and finally decays 
into hadrons. A schematic description of the pro- 
cess in light-cone kinematics is shown in FigfT] 
Boost-invariance can be justified in the central 
region of the collision since the observed distribu- 
tion of particles is flat, in agreement with the pre- 
diction of hydrodynamic boost-invariance, where 
(space-time) fluid and (energy-momentum) parti- 
cle rapidities are proven to be equal [5] (see later 
the discussion in section 3). 

The Bjorken flow was instrumental for deriving 
and predicting many qualitative and even semi- 
quantitative features of the quark-gluon plasma 
formation in heavy-ion reactions. It was later 
confirmed and developed through numerical sim- 
ulations where various effects have been included, 
up to the full 4-dimensional structure of the hy- 
drodynamic flow. However, since hydrodynamics 
is an effective mean-field approach, it says only 
little on the relation with the microscopic gauge 
field theory, i.e., in the present case, Quantum 
Chromodynamics (QCD). Some important ques- 
tions remain unsolved, such as the reason why the 
fluid behaves like a (almost) perfect fluid, what is 
the small amount of viscosity it may require, why 
and how fast thermalization proceeds, etc... 



The problem is made even more difficult by the 
strong coupling regime of QCD, which is very 
probably required, since a perturbative descrip- 
tion leads in general to high rj/s ratio. In fact, 
large collision cross-sections with transversally 
cut-off momenta are required in order to explain a 
small macroscopic shear viscosity. The mean free 
path and the average transverse momentum in- 
duced by the gauge field medium should be small, 
in order to damp the nearby force transversal to 
the flow measuring shear viscosity. In a classical 
physics context, one may write 
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where p par t, {resp. p a v) is the particle density, 
(resp. average transverse momentum) and A/ ree 
the mean free-path. Both small mean free-path 
and low average transverse momentum are ex- 
pected in the non-perturbativc regime of QCD 
and more generally from a gauge theory at strong 
coupling. We also quote in formula (Q]) the 
AdS/CFT value 1/4tt in terms of the funda- 
mental constants' ratio h/ks, corresponding to 
the large coupling limit of the conformal field 
theory. This AdS/CFT ratio appears to be very 
small compared to all known non-relativistic flu- 
ids, as suggested by the inequality between square 
brackets. Indeed, the AdS/CFT value has been 
proposed as an absolute theoretical lower bound 
[10] . A discussion is in progress to fully con- 
firm this statement or to find counter-examples 
in microscopic theories (see [11) and references 
therein) . 

In any case, a very small phenomenological vis- 
cosity ratio is confirmed; not meaning that vis- 
cosity effects should be neglected . Observations 
on the elliptic flow primarily, and other observ- 
ables, allow one to give estimates of rj/s which 
are within range of the AdS/CFT value (see, e.g. 
M). 

It is thus interesting to use our modern (while 
still in progress) knowledge of non-perturbative 
methods in quantum field theory to fill the gap 
between the macroscopic and microscopic de- 
scriptions of the quark-gluon plasma produced in 
heavy- ion collisions. Lattice gauge theory meth- 
ods are very useful to analyze the static proper- 
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ties of the quark-gluon plasma, but they are still 
unable to describe the plasma in collision. 

Hence we are led to rely upon the new tools of- 
fered by the Gauge/Gravity correspondence and 
in particular on the most studied and well-known 
of its realizations, namely the AdS/CFT duality 
[12] between N — 4 supcrsymmetric Yang-Mills 
theory and type IIB superstring in the large N c 
approximation. The properties of the gauge the- 
ory in (physical) Minkowski space in 3+1 dimen- 
sions arc in one-to-one relation with properties 
of the bulk theory, in the 10D target space of 
strings. For our applications, strong coupling fea- 
tures will be related by duality to the geomet- 
ric structure of the 5D metric with Minkowski 
boundary, which is built on top of vacuum AdS. 
At weak coupling, the physical space corresponds 
to the world-volume of a stack of infinitely many 
(N c — >oo) joint D3-branes. Note that the overall 
consistency of the correspondence scheme will be 
required, even if only geometrical properties will 
be used in practice. 

One should be aware, when using AdS/CFT 
tools, that there does not yet exist a gravity 
dual construction for QCD. Also, the relation 
between pcrturbative (weak-coupling) and non- 
perturbative (strong-coupling) aspects of QCD 
are left over in the correspondence, which con- 
veniently describes the strong coupling regime of 
gauge theories uniquely. However, the nice fea- 
ture of the quark-gluon plasma is that it is a de- 
confined phase of QCD, characterized by collec- 
tive degrees of freedom. Thus one may expect to 
get useful information from AdS/CFT. This has 
already been proved in the description of static 
geometries, by the evaluation of n/s [S]. Other 
aspects of the QGP can be studied via the dual- 
ity, in a static plasma configuration at fixed tem- 
perature. 

Our specific aim and the subject of the present 
lectures is the investigation of the Gauge/Gravity 
duality in a dynamical time- dependent setting us- 
ing the AdS/CFT correspondence. The goal is to 
describe a collision process in a strongly coupled 
gauge theory, including in particular the cool- 
ing of temperature accompanying the proper-time 
evolution of the fluid. 

In section [2 we introduce the theoretical tool 



of holographic renormalization and apply it, as 
a warm-up exercise, to a static uniform plasma. 
Sections [3] and |4] present, respectively, the late 
and early proper-time analysis of the plasma 
evolution, with particular emphasis on boost- 
invariant flows. 

2. AdS / CFT and strong gauge interactions 

In the previous section, we mentioned the ubiq- 
uity of hydrodynamic methods in the description 
of QGP produced at RHIC. Yet, despite their suc- 
cess in describing the data, we have to keep in 
mind that they are used as a phenomenological 
model, without a real derivation from gauge the- 
ory. This is understandable, since almost perfect 
fluid hydrodynamics is intrinsically a strong cou- 
pling phenomenon, for which one lacks a purely 
gauge theoretical methodic 

On the other hand, there exists a wide class 
of gauge theories, which can be studied analyti- 
cally at strong coupling. These are superconfor- 
mal field theories with gravity duals. In the unify- 
ing context of string theory, using the AdS/CFT 
correspondence, one is able to map gauge the- 
ory dynamics (CFT) at strong coupling and large 
number of colors into solving Einstein's equations 
in asymptotically anti-de Sitter space (AdS). The 
theories with gravity duals can differ substantially 
from real world QCD at zero temperature. The 
best known example of such theories - J\f = 4 su- 
per Yang-Mills (SYM) - is a superconformal field 
theory with matter in the adjoint representation 
of the gauge group SU(N C ). Because of the con- 
formal symmetry at the quantum level this theory 
does not exhibit confinement. On the other hand, 
differences between Af — 4 SYM and QCD are 
less significant above QCD's critical temperature, 
when quarks and gluons are in the deconfincd 
phase. Moreover, it was observed on the lattice 
that QCD exhibits a quasi-conformal trend [14] 
in a range of temperatures T > 300 MeV. There 
the equation of state is reasonably approximated 
by the conformal relation e = 3p, correspond- 
ing to a traceless energy-momentum tensor. The 

5 Lattice QCD methods do not work well here, as this 
would require analytical continuation to Minkowski signa- 
ture which is nontrivial in this context 1131 . 
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above observations, together with experimental 
results suggesting that QGP is a strongly coupled 
medium, are an incentive to use the AdS/CFT 
correspondence as a tool to get insight into non- 
perturbative dynamics. 

2.1. AdS/CFT correspondence 

We will now describe how to set up an 
AdS/CFT computation for determining the 
space-time behavior of the gauge field plasma 
corresponding to a Bjorken boost-invariant flow. 
The problem is formulated in terms of the proper- 
time evolution of the matter energy-momentum 
tensor on the CFT boundary [TSJ. This method 
docs not make any underlying assumption about 
local equilibrium or hydrodynamical behavior. 
We will indeed obtain the hydrodynamic expan- 
sion as a generic consequence of the late time 
behavior of the expanding and strongly coupled 
plasma. 

Suppose we consider some macroscopic state of 
the plasma characterized by a space-time profile 
of the energy-momentum tensor 



T^{x p ) , pt,v,p= {0,-- - ,3} 



(2) 



Then, since the AdS / CFT correspondence asserts 
an exact equivalence between gauge and string 
theories, such a state should have its counterpart 
on the gravity side of the correspondence. Typi- 
cally, it will be given by a modification of the ge- 
ometry of the original AdS§ x S 5 metric. This fol- 
lows from the fact that operators in the gauge the- 
ory correspond to fields in supergravity (or string 
theory) . When we consider a state with a nonzero 
expectation value of an operator, the dual grav- 
ity background will have the corresponding field 
modified from its 'vacuum' value. In the case of 
the energy-momentum tensor, the corresponding 
field is just the 5D metric. One then has to as- 
sume that the geometry is well defined, i.e. it is 
free from naked singularities - singularities not 
hidden by an event horizon -. This principle se- 
lects the allowed physical space-time profiles of 
the gauge theory energy-momentum tensor. Thus 
this becomes the main dynamical mechanism for 
selecting the appropriate strongly coupled gauge 
theory from the dual Einstein equations. In prac- 
tice, the geometry will be obtained from a "holo- 



graphic rcnormalization" procedure starting from 
the boundary. The boundary conditions serve as 
initial conditions for the construction of the bulk 
features, the fifth dimension acting as a renor- 
malization scale. Let us describe in detail this 
construction. 

2.2. Holographic renormalization 

The Gauge/ Gravity duality can be described 
as an "holographic" correspondence between the 
4-dimensional physical space where the gauge 
theory lives and the 5-dimensional space where 
the supergravity (weak curvature) approximation 
of the 10-dimcnsional string theory is valid. It 
means qualitatively that the whole information 
should be the same on both sides of the corre- 
spondence, despite the difference in dimensional- 
ity. In practice, this notion of "holography" has 
a precise realization in terms of the "holographic 
renormalization" program [16] . Let us illustrate 
the holographic rcnormalization by simple exam- 
ples, which will in fact be sufficient for the appli- 
cations we have in mind. 

Suppose that, when considering the presence of 
matter on the 4-dimensional physical space, our 
corresponding 5-dimensional geometry is param- 
eterized by 



j 2 g^(x p ,z)dx p dx v + dz 2 _ 



= gl D p dx a dx^ (3) 



where we adopt the Fefferman- Graham definition 
[T7] of the 5D metric. The flat case = rj^u 
parametrizes AdSs in Poincarc coordinates. The 
conformal boundary of space-time is at z = 0. 

Considering the general metric ((31) from the 
point-of-view of the AdS/CFT correspondence, 
the following questions are in order: 

i) What are the constraints imposed on 

9p.u{x p , z)l 

ii) What is the corresponding energy- 
momentum profile {T^ v {x p ))l 

The metric ([3]) has to be a solution of 5D Ein- 
stein's equation with negative cosmological con- 
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The expectation value of the energy momentum 
tensor may be recovered by expanding the metric 
near the boundary z — 0. The "holographic renor- 
malization" procedure |16j . for a Minkowskian 
boundary metric, give^l 



g liv {x p , z) = i)^ + z A g i $(x p ) + 
and 

N 2 



(5) 



(6) 



The relations (|5|6j) can be used in two ways. 
Firstly, given a solution of Einstein's equations, 
we may read off the corresponding gauge theoret- 
ical energy-momentum tensor. Secondly, given a 
traceless and conserved energy-momentum pro- 
file, one may integrate Einstein's equations into 
the bulk in order to obtain the dual geometrj@. 
Then the criterion of non singularity of the geom- 
etry will determine the allowed space-time evolu- 
tion of the plasma. Let us note that this formula- 
tion is in fact quite far from a conventional initial 
value problem. 

2.3. Example: static uniform plasma 

Before moving to the case of an expanding 
plasma, it is convenient to consider the simple 
situation of a static uniform plasma with a con- 
stant energy-momentum tensor. This will provide 
a useful exercise introducing the main tools for 
further applications. 



6 One can show that such solutions lift to 10D solutions 
of ten dimensional type IIB supcrgravity. The effective 
5D negative cosmological constant comes from the 5-form 
field in 10D supcrgravity. 

7 In principle, in J5j also additional logarithmic terms ap- 
pear (see |16p. They are absent in all the cases we are 
going to examine. 

8 This can be done order by order in z 2 , which is a near- 
boundary expansion. However potential singularities are 
hidden deep in the bulk, thus this power series needs to 
be resummed. 



Starting from a constant diagonal energy- 
momentum tensor 
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with pj| = p± and energy density and pressure re- 
lated by e = 3p, one has to solve the 5D Einstein's 
equations for an Ansatz 



-e a ^dt 2 + e b{z ^dx 2 + dz 2 



(8) 



obeying the above-mentioned boundary condi- 
tions (|5I6[) . Inserting (JSJ into the Einstein's 
equations Q and splitting them into compo- 
nents, one gets 

R n => -8a' + za' 2 - 66' + 3za'b' + 2za" = , 
i? 55 => -2a'+za' 2 -66' + 3z6' 2 + 2z(a"+36") = , 
i? 22 => -2a' - 126' + za'b' + 3zb' 2 + 2zb" = . 



Then 

-Rn — ^55 => a' 



zb 



/2 



2zb" 



zb 1 - 2 



(9) 



while combining it with the expression from R22 
for a' gives 



- 36'' 



zb" = . 



(10) 



From those last two equations, it is easy to ob- 
tain independent equations for a(z) and b(z) and 
their solutions. With appropriate normalizations 
at z = 0, one gets 



b EE l0g(l+-j 

z 



where zq is the integration constant. Hence, Ein- 
stein's equations can be exactly solved analyti- 
cally in this case and we find [15] that the exact 
dual geometry of such a system is 



ds 2 = - 



1- 



dt 2 
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This metric may look at first glance unfamiliar, 
but a change of coordinates 

(13) 



1 



transforms it to the standard metric form of the 
AdS-Schwarzschild static black holcH 



''ds 2 = - [ 1 - — ) dt 2 



dx 1 



1 



1 



-dz 2 , (14) 



with zq = zq/ s/2 being the location of the hori- 
zon. Before we proceed further, let us note here 
one crucial thing: the fact that the dual geome- 
try of a gauge theory system with constant energy 
density is a black hole was not an assumption, but 
rather an outcome of a computation. 
The Hawking temperature 



T 



1 

TTZ 



V2 
irz 



(15) 



is then identified with the gauge theory tem- 
perature, and the entropy with the Bekenstein- 
Hawking black hole entropy 



N 2 



_ N 2 T 3 
2 c 



(16) 



which is 3/4 of the entropy at zero coupling. This 
identification between the black hole characteris- 
tics and the equilibrium thermodynamics of the 
plasma reveals a striking correspondence between 
the near-horizon and boundary properties of the 
gauge/gravity system. 

To finish our discussion of the static black hole, 
we note that the Fefferman- Graham coordinates 
cover only the part of space-time lying outside 
of the horizon. It will appear useful to introduce 
different coordinates allowing one to go inside the 
horizon, such as the Eddington-Finkelstein's co- 



9 The identification of the metric (112 ft as the dual of a 
static fluid at fixed temperature has previously been made 
in 1181 using the approach 'from bulk to boundary'. This is 
the reverse way w.r.t. holographic renormalization which 
proceeds 'from boundary to bulk'. The interest of holo- 
graphic renormalization is that it allows a construction of 
the dual geometry starting only from 4-dimensional phys- 
ical data. 



ordinate system. Through the time coordinate 
redefinition 

1 _ / Z Zq + z\ 

u = t zq 2 arctan h log , (17) 

4 V z o zq-zJ 



the metric (fT4"]) becomes: 



~2 j 2 

z ds 



1 - — )du z 



dx 2 + 2dudz . (18) 



Now the metric is well defined at z — zq and 
the horizon can be smoothly crossed. Note that 
the time coordinate gets an infinite shift near the 
horizon, here being harmless, but which reveals 
to be quite subtle for the time-dependent back- 
grounds we will consider now. 

3. AdS/CFT and late time quark-gluon 
plasma flow 

In this part, we describe applications of the 
methods introduced in section[5]to the analysis of 
plasma dynamics in strongly coupled gauge the- 
ories. The time dependence of the system intro- 
duces considerable computational difficulties and 
it is in general a hard task to achieve a detailed 
quantitative description of the plasma evolution. 
Nevertheless, significant information on the dy- 
namics can be extracted working in simplified and 
symmetric, but realistic setups. We shall consider 
specific regimes of the problem in exam, but some 
general lessons can still be drawn about the highly 
nonlinear regime of the plasma expansion. 

Throughout the next sections we will model 
the products of heavy- ion collisions with a Af = 4 
SYM plasma at finite temperature. We mainly fo- 
cus on the results obtained in Ref. [15] concerning 
the late proper-time behavior of the quark-gluon 
plasma. Recent and detailed reviews on the sub- 
ject are also Refs. [19] and [20] . 

The structure of the sections is organized in the 
following way. In section [3Tl we review briefly the 
"Bjorkcn flow" description of the central rapid- 
ity region of heavy-ion collisions. This serves as 
an introduction to section 13.21 where we present 
the gauge theory setup which has been used in 
[15j to model the products of heavy-ion reac- 
tions. The description is in fact inspired by the 
hydrodynamic description of heavy-ion collisions 
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by Bjorken [8] and based on the central assump- 
tion of boost invariance. This hypothesis seems 
experimentally to be valid in the central rapidity 
region of the heavy-ion collision process. In sec- 
tion 13.31 we review the computation of the holo- 
graphic dual geometry of the plasma system in 
the asymptotic limit of large proper-time. The 
properties of this class of time-dependent space- 
times are investigated in section l3T4l The require- 
ment of the absence of naked singularities in the 
bulk selects uniquely the gravity dual of a perfect 
fluid. This non-singular solution can be inter- 
preted as a black hole moving off from the AdS 
boundary in the radial direction, as a function 
of proper-time. We conclude with sections 13.51 
and 13.01 where we present attempts of going be- 
yond the assumptions of perfect hydrodynamics 
and boost invariance. 



3.1. Hydrodynamic Bjorken flow 

The introduction of relativistic hydrodynam- 
ics in the description of high-energy hadronic re- 
actions is due to Landau in a remarkable pre- 
monitory paper [6]. However, the experimental 
evidence of a connection between heavy-ion col- 
lisions and relativistic hydrodynamics has been 
first modeled by Bjorken in [8]. In this descrip- 
tion, boost-invariance of the plasma arises as an 
outcome of the assumption of hydrodynamics, in 
a semi-classical approximation where the fluid fol- 
lows ballistic trajectories or, in other terms, the 
momentum-space rapidity is equal to the space- 
time one (see Figj2]). The model is therefore well 
suited to parametrize the central region of ultra- 
rclativistic collisions, for which the particle dis- 
tribution is nearly fiat as a function of rapidity. 

In this section, we review a number of aspects 
of this picture. The summary is intended to be 
useful in the follow-up to motivate the assump- 
tions of the AdS / CFT models we discuss and to 
interpret the physical results we obtain. 

For the discussion that follows, it is convenient 
to parametrize flat four-dimensional space-time 
of coordinates (t,xi,x±), where x± — (x2,xs), in 
terms of the light-cone coordinates (x + ,x~,x±). 
The latter are defined by x = t ± x\ , in which 
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Figure 2. Semi-classical fluid trajectories: the In- 
Out cascade. The rapidity y is equal to the space- 
time rapidity r\=\ log . 

the flat space-time metric reads 

ds 2 = -dx + dx" + dx\ . (19) 

The proper-time r and space-time rapidity rj of 
the fluid are related to the light-cone coordinates 
by x = t e ±v . 

In the first stages of the collision, a rather dense 
interacting medium is created and the individual 
partonic or hadronic degrees of freedom can be 
neglected to a good approximation. This obser- 
vation justifies the treatment of the medium as 
a fluid and, moreover, allows the assumption of 
local equilibrium. In the simplest perfect fluid ap- 
proximation, the energy-momentum tensor reads 

T» v = (e + p)u»u v -pr]^, (20) 

where e, p and are the energy density, pressure 
and 4- velocity of the fluid, respectively; while 
is the flat metric of expression ([T9| . We further 
assume an equation of state of the form 

e = gp, (21) 

where is the sound velocity of the liquid. In 

the case of a conformal theory, the traceless con- 
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dition Tjf = implies g = 3 in four dimensions. 
The continuity equation d^T^ = for a perfect 
fluid of energy density and pressure related by 
(l2~Tj) reduces to two equations 



gd+(\np) =- 



d+p- 



gd-Qnp) = ^^d^p-^e^d+y (22) 
for the pressure and the rapidity 



2 \e-p 



(23) 



of the fluid. The other thermodynamical quanti- 
ties can be derived from the equation of state (|2"Tj) 
and the standard thermodynamical identities 



p + e = Ts . 



de = Tds . 



(24) 



where T and s are the temperature and entropy 
density of the fluid and we have assumed, for sim- 
plicity, vanishing chemical potential. 
The result is 

p=-=p T 9+1 , s = s T 9 oc pi*T , (25) 
9 

for constant po and s . 

The key observation is that the expansion of 
the fluid can be described by a semi-classical pic- 
ture. In fact, the high occupation numbers of 
the system allow to assume that the fluid compo- 
nents follow quasi-classical trajectories in space- 
time, here straight-line trajectories starting at the 
origin (see FigJ5]), and given by 

y = r ] . (26) 

Plugging the Ansatz (|26| into (|22|) . one obtains 



gd+(\np) 



1 + 9 

2x+ 



<7<9_(lnp) 



from which it easily follows 

i P + - 

m — = log x x 

Po 2 



p(r) 



Pa 

1 + 8 

r a 



(28) 



The simple Ansatz (|2"6")l thus leads to a boost- 
invariant system. Such a physical description is 
suitable to model the central rapidity region of 
highly relativistic heavy-ion collisions, where a 



central plateau in the distribution of particles is 
detected. 

Specifically, in the perfect fluid four- 
dimensional case, the behavior is 

p(t) = ^ = ^—T 4 (t), S ~T 3 (t). (29) 

Note that the entropy density verifies s = sqt . 
Hence the total entropy S per unit of rapidity 
(and transverse area) is constant [8] (as expected 
from a perfect fluid) since 



dS — sd 3 x~sr / d 2 x±dy 



(30) 



3.2. Boost-invariant flow 

Inspired by Bjorken's analysis, the methods 
provided by the AdS/CFT correspondence have 
been applied to study the dynamics of boost- 
invariant strongly interacting gauge theory mat- 
ter in [15) . There are technical reasons for mak- 
ing this symmetry assumption, due to the com- 
plication of the Einstein's equations in the grav- 
ity dual. On a more theoretical level, using the 
M = 4 SYM theory as a substitute to QCD 
seems at first sight highly problematic. More- 
over, experimentally the assumption of boost in- 
variance is not optimal, since the observed mul- 
tiplicity distribution of particles in heavy-ion col- 
lisions is rapidity-flat only in the central region. 
In terms of the hydrodynamic modelization, this 
means that the entropy is rapidity- independent, 
and thus obeying boost-invariance, only in the 
same region. 

All these warning remarks have to be certainly 
taken into account, when using the AdS/CFT 
correspondence in a physical context. However, 
there are convincing arguments to show that for 
appropriate observables and by using AdS/CFT 
methods, one obtains a first fruitful approach to 
the behavior of QCD at strong coupling. From a 
phcnomcnological point-of-view, boost-invariance 
is widely used in modeling the relativistic heavy- 
ion collisions at RHIC in the central region of ra- 
pidity, where the hydrodynamic approach seems 
successful in describing the data and the occur- 
rence of the quark-gluon plasma. Moreover, this 
plasma is a non-confining phase of QCD for which 
no specific hadronic scale seems to be relevant. In 
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a more theoretical language, in this regime QCD 
appears to be nearly conformal. The above argu- 
ments suggest that it is reasonable to apply boost- 
invariance, at least in a first approximation, and 
to study the dynamics of the AT = 4 SYM plasma 
under consideration. 

In a boost-invariant setup, it is natural to 
parametrize the flat four-dimensional boundary 
in terms of the coordinates (r, rj,x±), where t is 
the proper-time, n the space-time rapidity and 
xj_ = (x2,X3) denotes the transverse directions. 
We already introduced these coordinates in the 
previous section, but it is worth being more ex- 
plicit here. In terms of the Minkowski coordinates 
(t, x\ 1 x±) 7 the definition of r and n reads 

f: = TC0sh77, X!=rsinh?7, (31) 

and the metric has the form 

ds 2 = -dr 2 + T 2 dif + dx\ . (32) 

For simplicity, further assume independence on 
the transverse coordinates x±, corresponding to 
the limit of infinitely large nuclei colliding. 

The energy-momentum tensor of the gauge 
theory can then be written in diagonal form and 
has only three non- vanishing components: T TTl 
T m and T X2X2 = T X3X3 , which depend on proper- 
time. The conservation condition £>„T MI/ = 
and, since we are dealing with a conformal theory, 
the tracelessness condition Tff = further restrict 
the form of the energy-momentum tensor. It can 
be shown that all components can be expressed 
in terms of a single function e = e(r): 

T^u =diag^,-r 2 (re' + e), e +^-,e+^-J . (33) 

The function e(r) can be interpreted as the en- 
ergy density of the plasma at mid-rapidity (i.e. 
at X\ = 0) as a function of proper-time. The last 
kinematic constraint on the stress-energy tensor 
is the positive energy condition. For any time- 
like vector the energy density in the reference 
frame whose time-like direction is specified by 
should be non-negative 

T„ u tn v > . (34) 

For the Ansatz (|33|). it implies 

e(r)>0, e'(r)<0, re'(r) > -4 e(r) . (35) 



The explicit expression of e specifies the dynam- 
ics of the theory and it is the information we aim 
to determine through the AdS/CFT correspon- 
dence. 

In the following, we restrict the analysis to the 
class of gauge theories whose late proper-time 
evolution is specified by a power-law parameter 
s, such that 

e(r) cx 1 . (36) 

The positive energy condition (|34l) fixes the range 
of possible powers s to < s < 4. 

Notice that the family of gauge theories de- 
scribed by ([55)1 . includes a number of cases of 
physical interest. If the longitudinal pressure T m 
vanishes, see (|33)l . the plasma is in a free stream- 
ing phase which is expected to be well suited to 
the case of a weakly coupled plasma [2T]. Solv- 
ing for e in (f3"3")) leads to e(r) cx r _1 . If the 
gauge theory plasma behaves, instead, as a per- 
fect fluid, the energy-momentum tensor has to 
be of the form (|2U)) with e = 3 p. Under the as- 
sumption of boost invariance, = (1, 0, 0, 0) and 
the comparison with (|3"3")) leads to the Bjorken 
solution of section I3TH e(r) cx r~ 4 / 3 . The case 
e(r) cx t° is also of interest, it describes a fully 
anisotropic medium with negative longitudinal 
pressure = T vv /t 2 = — e = —p±_ that might be 
relevant for the early times dynamics at strong 
coupling, as advocated in [22]. As we shall see in 
section [U the status of early times dynamics in 
the AdS/CFT correspondence seems to be more 
complex. 

While we will deal in detail with the dual de- 
scription of the early-time flow in the following, 
in the present part we restrict the analysis to the 
range < s < 4. In what follows, we will de- 
termine which, if any, of the above behaviors is 
relevant for the description of the large proper- 
time evolution. 

3.3. Large proper-time behavior of the 
gravity dual 

In this section, following the procedure of holo- 
graphic rcnormalization introduced in section [2~2l 
we repeat the construction of [15] of the dual ge- 
ometries to the boost invariant gauge theories 
with energy densities of the form (|3"6"|) . Since 
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the bulk metric shares the same symmetries of 
the plasma system, the most general Ansatz in 
Fcffcrman- Graham coordinates reads 

z 2 ds 2 = -e^d^+^e^dn 2 

+ e c{T - z) dx\ + dz 2 . (37) 

According to the holographic dictionary, to deter- 
mine the metric components cl(t,z), b(r,z) and 
c(t, z), Recalling Einstein equations with a nega- 
tive cosmological constant 

R*p = ^g a0 (R- 12) (38) 

one has to solve them in a power series expansion 
in the radial coordinate z. The indices a, j3 in (|38|) 
denote the bulk coordinates (r, z, ij, x±) and the 
metric components are subject to the boundary 
conditions 

o(r, z) = -z 4 e(T) + z 6 a 6 (r) + z 8 a 8 (r) + . . . (39) 

around z = and similar ones for b(r,z) and 
c(t,z). By studying the explicit form of the 
Fcffcrman- Graham expansion, it can be shown 
that the large r asymptotics of the exact solu- 
tions can be obtained analytically (see [15] for the 
details of the computation). In fact, introducing 
the scaling variable 

v = Ar, (40) 
r 4 

the exact solutions have an expansion of the form 

(41) 



a(T,z) = a(v) + o(Jz) > 



for some positive power a. The asymptotic equa- 
tions, obtained from (|3"B")) taking the limit r — > oo, 
while keeping v fixed, can then be solved exactly. 
The solution is 



a(v) — A(v) — 2m(v) 
b(v) = A(v) + (2s - 2)m{v) 
c(v) = A(v) + (2 - s)m(v) 
where 

A 2 (s)v s ) 
1 + A{s)v 4 



(42) 



A(v) 
m(v) 

A(«) 



5>"(i 



1 



In 



4A(s) 1 - A{s)v 4 



3s 2 - 8s 



(43) 



3.4. Selection of the perfect fluid 

Let us now examine in more detail the geome- 
try we obtained. The dual asymptotic metric of 
the boost-invariant plasma has a potential singu- 
larity for v = A _1 / 4 (s), where the argument of 
the logarithm in (|43[) vanishes. 

As usual, to distinguish a physical singularity 
from a coordinate one, it is sufficient to construct 
a scalar invariant formed out of the Riemann 
tensor which diverges at the singularity location. 
While the curvature scalar is everywhere regular, 
the simplest non-trivial scalar that can be con- 
sidered is the square of the Riemann tensor (or 
Kretschmann scalar) 

m 1 = R a ^ 5 R aPlS . (44) 

In the asymptotic limit r — > oo with v fixed, the 
general structure of this expression becomes 

W,s)_ MO (l. 



m 2 = 



(45) 



(1- A 2 (s)v 8 , 

for a positive j3. The explicit expression of the nu- 
merator N(v, s) is quite complicated and can be 
found in |15j . For generic s, the above asymptotic 
expression for *H 2 diverges at v = A _1 / 4 (s) and 
the singularity is a physical naked one. It turns 
out that only for s = 4/3 the fourth order pole 
of the denominator gets cancelled by an identical 
contribution in the numerator leaving a bounded 
scalar invariant. Only the physical singularity in 
the dual geometry of a perfect fluid (which will 
be found at z — oo, see later on) can be hidden 
behind an event horizon. 

The above analysis suggests a criterion to se- 
lect the profile e(r) in (|33|) . The proposal of [15] is 
that a physical boundary energy-momentum ten- 
sor has to be dual to a geometry satisfying the 
cosmic censorship. It then turns out that the re- 
quirement of the absence of naked singularities in 
the bulk uniquely selects perfect hydrodynamics 
as the large proper-time dynamics of the strongly 
coupled plasma. 

The scaling variable (|4"0"]) and the energy den- 
sity for a perfect relativistic fluid are 



— i e( T ) = — , 

T3 T3 



(46) 



24 



where we have reinstated a dimensionful param- 
eter eo in e. In this case, the metric coefficients 
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simplify considerably and the asymptotic geom- 
etry can be written in the compact Fcffcrman- 
Graham form 



W = - 



1 



LSI 

3 



1+ ^a- 

1 ^ 3 - 



dr 2 + dz 2 



(47) 



1 



e 7T 



(r 2 dy 2 



dx i 



3 r 4 /3 

This geometry is analogous to the static AdS 
black brane solution of section 12.21 and thus the 
potential singularity is indeed a coordinate one. 
Here, however, the horizon position depends on 
proper-time 



T3 



(48) 



Although a precise notion of temperature and en- 
tropy still lacks in a dynamical setup, in order to 
obtain a qualitative estimate one could naively 
generalize the static formula to infer 

r _ 4 (49) 

for the temperature and the total entropy S oc 
st = const., see (|30|) . Observe that these esti- 
mates are in agreement with the "Bjorken flow" 
expressions of equation (|29[) . Since the tempera- 
ture of the black hole coincides with the one of the 
plasma, the gravity dual of the plasma cooling off 
during expansion is in terms of a black hole which 
moves off in the radial direction of the bulk, from 
the boundary towards the interior. A similar pic- 
ture had already been qualitatively proposed in 

3.5. In-fiow viscosity and relaxation time 

It should be stressed that the results of section 
13.41 were obtained in the scaling limit. The grav- 
ity dual selects perfect hydrodynamics only in the 
strict asymptotic limit of infinite proper-time. In- 
deed, as shown for instance in [19] and [20], it is 
straightforward to repeat the computation of the 
square of the Ricmann tensor up to the first sub- 
leading corrections in the metric and check that 
it is always singular for perfect fluid dynamics. 
The evaluation of 9\ 2 leads to 

<R 2 = R (v) + \r 2 (v) , (50) 

T3 



where Rq(v) is finite, but R,2(v) develops a 4 th 
order singularity. From the gravity side, it is 
then clear that the energy density in (|46]) needs to 
be corrected. It turns out that these corrections 
are due to viscous hydrodynamics, but it is in- 
structive to see how the result arises using purely 
AdS/CFT methods. Assume that 



e(r) 



1 



1 - 



2% 



(51) 



for a generic positive exponent r. Solving Ein- 
stein's equations with this boundary condition 
and computing the Kretschmann scalar yields a 
result of the form 



D\ 2 = R Q (v) + \Ri(v) 



\Ri{v)i 

T3 



(52) 



where the last two terms are always singular. To 
obtain a bounded Kretschmann scalar, the only 
possibility is to make those two terms propor- 
tional to each other and therefore to set r = 2/3. 
This is exactly the scaling of a viscosity correc- 
tion to Bjorken hydrodynamics with shear viscos- 
ity 7/ = t]o/t (which follows from n oc T 3 in the 
N = 4 SYM case). Moreover, in PI] it has been 
shown that, to implement the cancellation of the 
potentially singular terms, the coefficient r]Q has 
to be tuned to 770 = 2" 1 / 2 3" 3/4 . This specific 
value corresponds to a shear viscosity to entropy 
ratio 



1 

47T 



(53) 



It is remarkable that the above result, obtained in 
a non-linear and dynamical setup coincides with 
the shear viscosity value obtained in [S] , studying 
the response of a static plasma at fixed tempera- 
ture to small perturbations. 

The procedure can be further generalized to 
higher orders to compute the other hydrodynam- 
ics coefficients. In general, the metric components 
a(r, z), b(r, z) and c(r, z) in (|3T|) have an expan- 
sion of the form 



1 



(54) 



12 



A. Bcrnamonti and R. Pcschanski 



while the one of the curvature invariant is 
$R 2 = ^ii„-L- (55) 



T 3 



The 2 nd order formalism is obtained truncating 
the sums above at n = 3 and has been worked 
out in [55] (with a correction given in Ref. |26j . 
see further on). It results into a relaxation time 



2 - In 2 
2ttT 



(56) 



which is about thirty times shorter than the one 
estimated on the basis of Boltzmann's kinetic the- 
ory. The value ([56]) has been corrected by match- 
ing the results for the Bjorken flow of Ref. [25] 
with the full second-order formalism for hydrody- 
namics in conformal field theories at finite tem- 
perature derived in [26] . 

Surprisingly, at this order, a new feature arises, 
namely the presence of a leftover logarithmic sin- 
gularity in the expansion of D\ 2 . In [25J, it has 
been argued that such a divergence could be cured 
considering the coupling of a dilaton to the bulk 
metric. A non vanishing dilaton profile turns on, 
on the gauge theory side, a non-zero expectation 
value of TrF 2 , meaning that electric and mag- 
netic modes are no longer equilibrated. The tun- 
ing of the dilaton field to achieve a non-singular 
Krctschmann scalar, leads to a negative expec- 
tation value (Tr_F 2 ) < 0, implying that electric 
modes dominate. 

However, the logarithmic divergences that arise 
at NNLO in the energy density expansion, have 
been considered more extensively in [57] and [55J 
(see also [55] and [3D])- The conclusion is that the 
logarithmic singularity origins in a pathology of 
the Fcfferman- Graham coordinates. There exists 
in fact a singular transformation to Eddington- 
Finkelstein coordinates, which yields a regular 
and smooth metric from the boundary up to the 
standard black brane singularity. In this coor- 
dinate system, the late proper-time expansions 
of curvature invariants is regular at the horizon. 
In order to proceed to higher orders, one would 
then have to perform the analysis in Eddington- 
Finkelstein coordinates. 



3.6. Duality beyond boost-invariance 

The efficacy of the approach that we presented 
in the previous sections resides in that no assump- 
tion about the dynamics of the plasma has to be 
made. It is the gravity dual analysis alone that 
selects the physically relevant evolution. How- 
ever, in order to solve the full system of Einstein's 
equations one is in general constrained to work in 
highly symmetric setups and in specific regimes of 
interest. In particular, it is very hard to relax the 
assumption of boost invariance of the gauge the- 
ory matter. Significant progress in this direction 
has been achieved through a different approach 
in [31] . where it has been shown in general how 
hydrodynamics arises from the gravity side. 

The starting point of the analysis is the five- 
dimensional static AdS black brane which is dual 
to a strongly coupled plasma at rest (i.e. with flow 
vector = (1, 0, 0, 0)) and at temperature T. In 
incoming Eddington-Finkelstein coordinates, the 
metric reads 



ds 2 



1 - 



dV + 2 dvdr 



+r 2 r]ij dx l dx l , 



(57) 



where the coordinate r is related to the radial 
coordinate z by r = 1/z. The conformal bound- 
ary of space-time is at r = oo and the horizon 
at r = T/ir. This choice of the coordinate sys- 
tem is consistent with the comments of the end 
of section 13. 5> these coordinates are well defined 
on the horizon and extend all the way from the 
boundary to the black brane singularity. 

By performing a boost, one can obtain the dual 
geometry to a uniformly moving plasma with 4- 
velocity m m 



ds 2 



i ^ 

7rr 



u^Uy dx^dx v (58) 



-2u^ dx^dr + r 2 (n^ u + u^u u )dx^dx 1 ' 



The main idea of [31j is to promote T and u 11 
to be slowly varying functions of the boundary 
Minkowski coordinates. In this way, the geometry 
ceases to be an exact solution of Einstein's equa- 
tions because of the appearance of gradients of 
the temperature and of the velocity. The metric 
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can be corrected order by order in a deriva- 
tive expansion. The integration constants that 
arise at each order can be fixed requiring regu- 
larity of the metric at the horizon. The energy- 
momentum tensor of the plasma can then be read 
off from the corrected geometry and it is ex- 
pressed, as the metric itself, in terms of T, 
and their derivatives. Explicit results have been 
worked out up to second order and read 



1 



16ttG 5 



(7rT) 4 (^ +4u"u") 



(59) 



-2(ttT) V tl/ + (ttT^) [ log 2 



-(2-log2) l-T£ 



rpliV 

~- L 2d ' 



where G5 is the five dimensional Newton's con- 
stant and the rather lengthy expressions of the 
quantities o^" , T%£, etc. can be found in the orig- 
inal paper [31] . The first term in T^ v is simply the 
perfect fluid energy-momentum tensor, the sec- 
ond term involves the shear viscosity, while the 
third is the contribution of second order hydro- 
dynamics. 

The conservation equation d^T^ = is equiv- 
alent, by definition, to the relativistic Navier- 
Stokes equation. Therefore, the dual geometry 
of every solution of the hydrodynamic equations 
is given by ([58]) plus correction terms up to the 
same order in the gradient expansion. 

We should note, however, that this construc- 
tion is adequate exclusively for near the equilib- 
rium setups. Starting from the boosted black 
brane metric, in fact, one implicitly assumes 
the existence of a hydrodynamical description, in 
terms of flow velocity and energy density. As we 
explain in section 0] when dealing with processes 
that do not admit such a description, as the early 
stages of heavy-ion collisions, one has to recur to 
different methods. 

4. AdS/CFT and early time quark-gluon 
plasma flow 

One of the main theoretical challenges of the 
physics of quark-gluon plasma is to describe its 



thermalization process. 

The plasma is created in the initial stages of 
heavy-ion collisions, during which the system is 
certainly out-of-equilibrium. However, simula- 
tions based on the RHIC data suggest that an 
hydrodynamical description becomes reliable at 
times < 1 fm/c after the collision event |32|33|34] . 
Pcrturbativc QCD calculations [33135136] lead to 
a much longer thermalization time Ttherm ^2.5 
fm/c, pointing in the direction of a strongly cou- 
pled scenario. As we saw, holographic methods 
are, up to now, one of the few tools of investi- 
gation available in this context, as well as a very 
promising one. 

In this section, we describe some steps for- 
wards towards a description of the thermaliza- 
tion process of QGP, modeled by a strongly cou- 
pled Af = 4 SYM plasma at finite temperature. 
Since it is a developing subject on a settled theo- 
retical question dealing with far-from-equilibrium 
dynamics at strong coupling, we will discuss three 
different approaches. 

i) In section 14.11 we presents estimates of the 
thermalization response-time computed through 
the quasi-normal modes of small perturbations 
away from the thermal background. This explores 
the information one can get on thermalization 
and isotropization from linearized solutions, cou- 
pling non-hydrodynamic perturbations to a mov- 
ing black hole background. 

ii) Section 14.21 is devoted to the discussion 
of the fully non-linear thermalization dynamics 
of a boost-invariant plasma, evolving from an 
anisotropic out-of-equilibrium early stage towards 
the hydrodynamic regime. The results point 
towards a dependence on initial conditions but 
also to the indication of a fast (if not complete) 
isotropization as characteristic of strongly cou- 
pled dynamics. 

iii) In section 14.31 we mention the results ob- 
tained in modeling the Lorentz-contracted rela- 
tivistic nuclei of a heavy-ion collision through ini- 
tial shock waves and look for a solution of the dual 
gravitational field after the collision. This pro- 
vides a realistic scenario, at least from the kine- 
matic point of view. 

A warning is in order to remind the follow- 
ing unsolved problem. The preferred phenomeno- 
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logical scenario for the QGP formation starts 
from initial conditions which may be dominated 
by weak QCD coupling, due to the high den- 
sity of partons in a fast boosted nucleus. Hence 
the problem remains open to describe on a mi- 
croscopic basis the transition from a far-from- 
cquilibrium weakly coupled system towards a hy- 
drodynamic one at strong coupling. At least the 
AdS/CFT correspondence allows to make a few 
first (and quite unique at present) steps on the 
understanding of strongly coupled dynamics. 

4.1. Thermalization response-time 
4.1.1. Quasi- normal modes 

The quasi-normal modes of a black hole de- 
scribe the response of the system after a small 
perturbation. Since the field that sources the per- 
turbation can fall into the black hole or radiate 
to infinity, the modes of oscillations decay, their 
frequencies being complex. The geometry there- 
fore undergoes damped oscillations that produce 
an exponential decay of the perturbation. The 
frequencies and damping of these oscillations are 
entirely fixed by the black hole, and are indepen- 
dent of the initial perturbation. For black holes in 
asymptotically flat space-times, the exponential 
decay of the quasi-normal modes is followed by 
a power-law decrease [37]. The decay is instead 
purely exponential in the case of black holes in 
AdS [35]. 

In the AdS/CFT framework, perturbing a 
static black hole in AdS corresponds to perturb- 
ing the approximately thermal dual state, and the 
decay of the perturbations describes relaxation 
back to thermal equilibrium. Computing the 
quasi-normal modes of the equations of motion 
linearized around the background allows there- 
fore to estimate the thermalization response-time 
of the strongly coupled dual gauge theory after a 
small deviation from equilibrium. 

The behavior of quasi-normal modes should 
be contrasted with the one of hydrodynamical 
excitations that go to zero at small transverse 
momentum and have a slower power-law de- 
cay. In the static case, for instance, both quasi- 
normal modes and viscosity calculation corre- 
spond to poles of specific retarded propagators, 
respectively for a scalar and metric deformation, 
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through Fourier transform in the time variable. 
However, in the limit of zero transverse momen- 
tum, the poles corresponding to the viscosity case 
vanish. The "non-hydrodynamic" nature of the 
scalar quasi-normal modes allows instead to ob- 
tain a finite thermalization response-time in this 
limit. 

The quasi-normal modes analysis has been car- 
ried out in [3J|] for black holes in four, five and 
seven dimensional global AdS and perturbations 
produced by a minimally coupled scalar field exci- 
tation. See also |40|41|42] for later developments. 
In asymptotically AdS space times, quasi-normal 
modes are defined as solutions of the wave equa- 
tion that are purely ingoing near the horizon and 
vanish at the boundary. For a minimally coupled 
scalar perturbation, they are obtained solving the 
wave equation for a masslcss scalar field V 2< J> = 
in the gravitational background, with absorbing 
boundary conditions at the horizon z = zq and 
Dirichlet conditions at the boundary z = 0. One 
writes 

V 2 $ = -^=d Q (V=gg af3 d p <S>) = (60) 

for the standard covariant coupling of the scalar 
field <3> to the background metric g a @ of determi- 
nant g. 

4.1.2. Static black hole 

Since the main example throughout the notes 
is AT = 4 SYM in four-dimensional Minkowski 
space, in the following we refer the results ob- 
tained in the case of a static planar black hole 
in AdS5. To facilitate the comparison with the 
results of [J3J discussed later in the section, we 
work in Feffcrman- Graham coordinates. 

Recall that in these coordinates, see (|T2]) . the 
metric of a five-dimensional black brane reads 

fi+gj z V z oJ z z 

where (t, x) are the boundary coordinates, z is 
the AdS radius and Zq is the location of the event 
horizon in the bulk. 

In the case of the static metric (joTj) . Eq. foT)]) 
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writes 



1 - 



z 3 1 



d 2 Mt,z) 



3 (1-3) ».*(*.*) 



. (62) 



(63) 



the wave equation (|62j) leads to a Heun equation 
in the bulk variable |40] 



+ <9* 

Using the separation of variables 
*(t t z) = e iut <f>(z), 



1 - S 



:2 



z(l-z)(2-z) 



(ZOLJ)' 



85(1 - z){2 - 2) 



(64) 



where z = (1 - (z/z ) 2 ) 2 /(l + (z/z ) 4 ), z = 
V^ttT)" 1 , cf. (JTSJ, and the prime denotes the 
derivative with respect to z. 

The quasi-normal modes analysis then pro- 
ceeds through a combination of analytic and nu- 
merical methods (see [39] for details). Note that 
quasi-normal modes appear also as poles in the 
complex frequency plane of the Fourier transform 
in time t of the retarded Green's function (see e.g. 

W)- 

Interestingly enough, for a homogeneous scalar 
excitation, the dominant decay mode at large 
time is given by a quasi-normal mode whose fre- 
quency acquires a non zero imaginary part and 
thus corresponds to a minimal exponential decay 
mode with temperature. Technically, by match- 
ing Eq. (|64p with a Shrodingcr equation 



dl,6(z*) + 



V8e 12z * /uj_y 



sinh 3/2 (8z*) WlV 
through the change of variable 



f ) = (65) 



z^z" = - [tanh-V/*o)] 4 . 

one gets quantized "energy states". The dom- 
inant exponential decay mode, corresponds to 
the state with minimal energy in the Shrodinger 



equation (|65|) . Denoting u> k the eigenvalues of Eq. 
(|65p . one finds 



— = a k -ib k , 
and 



>b k > 



>h = 2.74667, 



(66) 



(67) 



where T denotes both the plasma and black 
branc temperature. The corresponding response 
timescale is thus given by the first eigenvalue oj\ 
and leads to T resp = 1/Smwi ~ 0.116/T. 

4.1.3. Moving black hole 

A similar analysis was carried out in [43] for 
a black hole moving away in the radial bulk di- 
rection. This geometry was discussed in section 
13.41 and is the gravity dual of a boost-invariant 
expanding perfect fluid in Af — 4 SYM at large 
proper-time. Knowledge of the quasi-normal 
modes of this black hole gives an estimate of the 
thermalization response-time of the strongly cou- 
pled gauge theory, that is the relaxation time af- 
ter a perturbation of the local thermal equilib- 
rium due to the coupling to a scalar field. 

The background geometry describing a black 
hole moving off from the boundary in the z di- 
rection was given in (|48|) . In terms of the scaling 
variable v = zt" 1 / 3 , it writes 

2,2 (l-« 4 ) 2 j2 

z 2 ds z = --, irdT 1 + 

(1 + v 4 ) 

+ {l + v i ){r 2 dy 2 + dx 2 L ) + dz 2 . (68) 

Consider now the wave equation of a canonically 
coupled scalar field $ in the boost-invariant set- 
ting. In the scaling limit v =const and t — > 00, 
the corresponding Klein-Gordon equation (|60|) is 
now expressed in terms of the two variables r 
and z. Through the change of variables (r, z) — > 
(r, v), so that 



d z -> t *d v , d T -> d T - -t 3 d v , 



(69) 



and neglecting the non diagonal terms in the large 
r expansion, one gets 

1 (l+« 4 ) 2 „ 2 ,,, . 



+ t 3 d v 



;(l-v s )d v <£(T,v) =0 , (70) 
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which is similar to (|62[) but with a r-dcpcndcnt 
coefficient. Performing the separation of variables 



$(t,z) = /(t)0( W ), 

it reduces to two decoupled equations 



■ UJ 2 T 



(1-v 8 ) 



d v 4>{v) 



= 

, 2 (1-^ 4 ) 2 

v 3 (l+v 4 ) 



(71) 

(72) 
= 0. 
(73) 



where we introduced an a priori arbitrary param- 
eter i2>, whose allowed values will be fixed by the 
(quantized) solutions of (|73|) . 

Equation (|72|) is solved by linear combinations 
of the Bessel functions 



VtJ+3 



-wt 3 



with relevant large r behavior 



f(r) 



r 6 e 2 



i([ii)r3 



(74) 



(75) 



xsm uji = 2.74667 is the same as in (|67p . 
Indeed, the equation in the v variable ((75)) is for- 
mally identical to the Heun equation (|64p of the 
static black hole and therefore leads to the same 
numerical values for the quasi-normal modes and 
frequencies. However, the variables v and r of 
((71")) arc different from the variables z and t rele- 
vant in the static case. Moreover, the proper-time 
dependence (|75[) has a non trivial scaling com- 
pared to the plane wave dependence of the static 
solution. 

Noting that the moving horizon is at z = 

i 

4 t3 = ^p, cf. (|48l49p . it was suggested in 

(43) , that this result could be understood in terms 
of an adiabatic approximation of the expanding 
case, with locally fixed temperature. Indeed, 



l/(r)| 



|Sm(iji) t3 



-tt&iTt 



(76) 



where u)\ = 7rT(ai — ib\) is the dominant decay 
mode in the static case ((67)) . 

Plugging in the numerical value of the domi- 
nant decay mode (jB"7)) in ((75)) leads to a damping 
factor of the form 

exp ( - | • 2.7466 ■ r* ) ~ e~ 8 - 3TT . (77) 



The above results were extended in [41] to include 
also vector and tensor quasi-normal modes. 

In the limits to which this analysis can be con- 
sidered quantitatively relevant for the thermaliza- 
tion in heavy-ion reactions, the fast decay modes 
of frequency (|67p damp out to 1 / e of their original 
amplitude in a time no greater than 

i0.08fm/c, (78) 



T=T V 



T 



T=T„, 



for a typical initial peak gauge theory tempera- 



ture TL 



300 MeV, which is commonly con- 



sidered as reasonable in phenomenology. Given 
the highly anisotropic momentum space distri- 
bution expected in the early stages of a RHIC 
collision, one tentatively estimates a certain fi- 
nite number n of e-foldings of the "response ther- 
malization process" have to elapse before hydro- 
dynamic approximations can be used. The esti- 
mated thermalization time becomes 



Tthe 



n x t. 



resp 



T=T„, 



(79) 



leading to T t herm ~ 0.3 fm/c following the guess 
of [JT] of about n = 4 iterations. The above es- 
timate compares in order of magnitude with the 
simulations predictions on real experiments. Nev- 
ertheless, it should be clear that it is extrapolated 
from the quasi-normal modes analysis, that, by 
definition, only applies to the late-time stages of 
the whole thermalization process. Moreover, the 
thermalization response time has been obtained 
in a linearized approximation around the black 
hole background through AdS/CFT methods. 

As a final remark on this response time stud- 
ies, it is worth noting that the same properties 
are valid for quasi-normal modes dominating the 
transverse space diagonalization of the metric af- 
ter a non diagonal perturbation. Let us indeed 
introduce a non diagonal transverse component 
gx!,x 2 using its contravariant expression 

..2 



1 



(80) 



It verifies (with (z, r) = /(r) g%\ (v)) the same 
equations (|72|73j) : such as 



{l-v s )d v g*Uv: 



Gj 2 (l + i; 4 ) 



4\2 



3 1 4 

,jA I _ yl 



(81) 
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In some sense, isotropization and thermaliza- 
tion response-times are equal at large proper- 
times. We now will turn to the study of the non- 
linear problem of isotropization and thcrmaliza- 
tion when starting from far-from-equilibrium ini- 
tial conditions. 

4.2. Thermalization of boost-invariant 
plasma 

In the present section we describe three differ- 
ent approaches to the study of the thermaliza- 
tion process of a strongly coupled boost-invariant 
plasma. The analysis are to be intended as com- 
plementary and point towards a fast isotropiza- 
tion/thcrmalization. However, as we shall dis- 
cuss, this represents only an early stage for a full 
understanding of far-from-equilibrium dynamics 
leading to the QGP formation. 

4.2.1. Solution with "full anisotropy" 

One of the first holographic estimates of the 
isotropization/thcrmalization time of the quark- 
gluon plasma appeared in [22] . The setup of this 
analysis is similar to the one considered in [IB] , 
but with differences which we will comment at 
the end of this subsection. 

The starting proposal of the paper is that 
Bjorken hydrodynamics at late times can be sin- 
gled out, in the holographic renormalization pro- 
gram, relaxing the assumptions to only require 
that the metric tensor is a real and single-valued 
function of the coordinates everywhere in the 
bulk, without imposing any constraint on the 
curvature invariants. Applying the same strat- 
egy to early time dynamics |22| leads to infer 
the existence of a solution corresponding to the 
fully anisotropic case with constant energy den- 
sity e. In the notation of section 13.21 it corre- 
sponds to the solution with s = in equation 
([36)) , and thus to "full anisotropy" at initial time, 
namely e = p± = — pii . The system is initially 
anisotropic, with negative longitudinal pressure, 
and is expected to evolve at late times to ideal 
Bjorken hydrodynamics, with longitudinal and 
transverse pressure components equal and posi- 
tive. Isotropization must take place at some in- 
termediate proper-time. 

Assuming a smooth transition between full 



anisotropy and full isotropy, the isotropization 
time r, so can be estimated matching the small 
and large proper-time regimes. In [22], this 
proper-time was defined as the crossing value of 
the branch-point singularities of the two regimes, 
leading to 



Here cq is the same dimensionfull parameter that 
appeared in (|46|) in the late proper-time behavior 
e(r) = e r- 4 / 3 . 

Phenomenologically, to obtain a rough estimate 
of the isotropization time (|82p . one can extrapo- 
late the realistic evaluations of RHIC to the su- 
persymmetric case under consideration. The hy- 
drodynamical simulations for central Au+Au col- 
lisions with y/s = 200 GeV yield the energy den- 
sity e = 15 GeV/fm 3 at the proper-time r = 0.6 
fm/c. Plugging in the corresponding value of eo 
in (|82p and setting AT = 3, in [22] it was obtained 
t iso ps 0.3 fm/c. 

Ref . [23] thus points towards a small value of the 
thermalization time. However, the analysis relies 
on the assumption of a rapid and smooth inter- 
polation between the initial fluid conditions and 
the late time regime. A full solution of that tran- 
sition remains necessary, as we shall discuss now, 
either by starting from given initial conditions or 
by a general study of the non-linear gravitational 
solution in the bulk of the AdS space. 

4.2.2. Far-from-equilibrium forcing dy- 
namics 

A way to prepare an out-of-equilibrium state 
is to turn on, in the ground state of a system, a 
time-dependent perturbation of the background 
fields. For instance, one can introduce a time- 
dependent deformation of the boundary geome- 
try with compact support. From an holographic 
perspective, this acts as a source of gravitational 
radiation in the bulk, leading to gravitational col- 
lapse and inevitably to black hole formation. Af- 
ter the perturbation is turned off, the geometry 
relaxes to a smooth and slowly varying form. In 
the dual theory, the latter process describes re- 
laxation of the non-hydrodynamic degrees of free- 
dom towards perfect fluid-dynamics. 
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This approach has been taken in [55] (follow- 
ing previous analysis |46j ) to study the thermal- 
ization process of a boost-invariant TV = 4 SYM 
plasma. 

The familiar Ansatz of two-dimensional spa- 
tial homogeneity, 0(2) rotational invariance in 
the transverse plane (defined by x±_) and boost- 
invariance in the longitudinal direction (x^ = x\) 
has been considered in [35]. This setup has the 
advantage of allowing the comparison with the 
late proper-time results of [15] . Since the late hy- 
drodynamic evolution is nearly isotropic, it is es- 
pecially interesting to generate a high anisotropy 
in the initial non-equilibrium state. 

To get some intuition on the process, we borrow 
Fig|3]from [45]. The figure depicts a space-time 
diagram of the field theory evolution. The sys- 
tem starts out at proper-time r = in the ground 
state. At t = Ti, the four-dimensional boundary 
geometry starts to deform in the vicinity of x» = 
and the perturbation acts until time Tf, propa- 
gating in the ixy directions at speeds asymptot- 
ically approaching the speed of light. The rele- 
vant space-time region is I in red. At r = Tf, the 
field theory is out-of-equilibrium and significantly 
anisotropic. It relaxes towards local equilibrium 
and hydrodynamics in region II in yellow. After 
time r* , region III in green, the system is well ap- 
proximated by a perfect fluid behavior. Deriving 
the thermalization time of the plasma amounts 
to the computation of the proper-time in the 
diagram. 

In [45], the deformation of the boundary ge- 
ometry was taken to be produced by a time- 
dependent shear 



ds 2 = -dr 2 
Here 

7(t) =ce[i(T)]i(T)V^, 
with 

(r-T )\ 



r 2 e- 2 ^d 



If 



5(r) = 1 - 



A 2 



(83) 



(84) 



(85) 



c is a constant characterizing the amplitude of the 
perturbation and is the unit step function. The 
deformation acts during the time interval (t^, Tf), 
with endpoints Ti = tq — A and Tf = tq + A. 



Ti 



Figure 3. Space-time diagram of the field the- 
ory plasma evolution. The metric deformation 
is turned on in the time interval (t,;,t/). At r* 
hydrodynamics sets in. Figure from |45j . 



Choosing tq = 5/4 A, the geometry is flat at t = 
0. All quantities will be measured in units where 
A = 1, so that n = 1/4 and r/ = 9/4. 

The most general bulk Ansatz consistent with 
diffcomorphism and spatial 3D translational in- 
variance, together with 0(2) rotation invariance, 
is 



ds 2 



-Adr 



IdrdT 



B dx 2 



-2B 



dy 2 



(86) 



where A, B and S are all functions of the bulk 
radial coordinates r and t only. The coordinate 
r is related to the previously used z, by r = 1/z. 
Notice that here r and r are generalized infalling 
Eddington-Finkelstein coordinates. 

Einstein's equations, subject to the boundary 
conditions that the boundary metric g^ v (x) coin- 
cides with (|83|) . can be solved numerically. The 
boundary stress-energy tensor can be computed 
through 



SS G 



y/-g(x) 



(87) 



where Sq denotes the gravitational action. The 
time-dependence of the energy density, longitu- 
dinal and transverse pressures then follows from 
(see [16] ) 



N 2 

T " = 2 _ ^ diag ( _e ' p H' p - L ' ?5 - L ) • 



(88) 
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T — Tf T — Tf T — Tf 



Figure 4. Energy density (blue top curve), trans- 
verse pressure (green middle line) and longitudi- 
nal pressure (red bottom), divided by N^/2tt 2 , 
as a function of time for c = 1/4, c = 1 and 
c = 3/2. The dashed black lines show the sec- 
ond order viscous hydrodynamic approximation 
to the different stress tensor components. 



The numerical solution explicitly shows the for- 
mation of an event and apparent horizon, which 
build up in the bulk and merge at late times. At 
asymptotically late times, the metric settles down 
to the moving black hole dual of the perfect fluid 

The transition from the far-from-equilibrium 
behavior to near-local-equilibrium hydrodynam- 
ics is found to be regulated by the relative 
importance of the exponentially decaying non- 
hydrodynamic degrees of freedom in comparison 
to the slowly relaxing hydrodynamic modes. This 
is in contrast with the naive expectation that hy- 
drodynamic should set in at the time at which 
higher order terms in the hydrodynamic expan- 
sion become comparable to lower order terms. 

A quantification of the results is provided in 
figure FigJH which was presented in [45]. The 
energy density (e), transverse and longitudinal 
pressures (p± and py) are plotted for different 
values of c, function of time, starting at 

t = Tf. On top of the numerical data are rep- 
resented the hydrodynamic approximations. The 
plots show a significant anisotropy, even at late 
times, where hydrodynamics applies. For all val- 
ues of c, the transverse pressure approaches the 
longitudinal pressure from above. This statement 
can be checked by perturbative computations of 



the first order viscous corrections. As the value 
of c increases, the system gets closer to equilib- 
rium at Tf. Indeed, for larger c, the perturbation 
docs more work on the geometry and the system 
reaches a higher temperature. In a conformal the- 
ory, as J\f — 4 SYM is, the relaxation time of the 
non-hydrodynamic degrees of freedom must scale 
inversely with the temperature. For c — > oo, the 
system will always be close to local equilibrium 
and isotropic already at Tf. 

The thermalization time was defined in 
[45] as the time beyond which the stress ten- 
sor agrees with the hydrodynamic approximation 
within 10%. In all the analyzed cases, the relevant 
dynamics, from the production of the plasma to 
his relaxation to near to local equilibrium, takes 
place over a time — r, < 2/T*, where de- 
notes the local temperature at the onset of the 
hydrodynamic regime. 

4.2.3. General approach to the early-time 
flow 

In this section we discuss the results obtained 
in [47] , where an analysis of the early time expan- 
sion of the boost-invariant flow has been carried 
out on the lines of [15]. We refer the reader back 
to section [3~2l for the setup and the basic assump- 
tions. 

The idea of the approach is the same as for 
the large proper-time regime discussed in sections 
13.31 and 13.41 In the following, we describe how to 
holographically reconstruct the dual geometry to 
a plasma configuration with stress-energy tensor 
(|33l) at early proper-times. This is accomplished 
solving Einstein's equations in a power series ex- 
pansion starting from the conformal boundary of 
space-time and subject to the appropriate bound- 
ary conditions. 

In the late times case, the structure of Ein- 
stein's equations naturally led to the introduction 
of a scaling variable and to a remarkable simpli- 
fication of the computations. It turns out that 
it is not possible in general to apply the same 
procedure to the small r regime. In what follows 
we schematically reproduce the main steps of the 
computation and point out the complications en- 
countered in [47] . 

We recall that the most general Ansatz for the 
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bulk metric consistent with the plasma symme- 
tries is 

z 2 ds 2 = 



-e< T ^dr 2 +T 2 e b(T ^dy 2 
+ e c{T ^dx 2 L + dz 2 . 



Solving Einstein's equations ([38 

R Q p = \g a e{R- 12) 



(89) 



(90) 



in a power series expansion in the bulk coordinate 
z, together with the boundary conditions (|39[) 

a(r, z) = -z 4 e(r) + z 6 a 6 (r) + z 8 a 8 (r) + . . . (91) 

allows to fix the metric components a(r, z), b(r, z) 
and c(r, z) order by order in the z expansion. 

Analyzing the structure of the early-time ex- 
pansion of the equations, in [47] it was concluded 
that a scaling variable does not exist in general 
in the r — > limit. As we shall see, the solution 
of [35] is not unique. One has to analyze the full 
solution at r ~ 0. 

From a physical point of view, it is natural to 
expect that the initial conditions should play a 
crucial role at early times. On the gauge the- 
ory side, the initial state can be prepared in a 
multitude of ways and therefore its early-time be- 
havior cannot be expected to be universal. This 
contrasts the late time dynamics which, under 
the specific symmetry assumptions under con- 
sideration, is governed by a single temperature 
scale. The dissipative effects which act during 
the plasma evolution wash out the differences due 
to the initial conditions leading to a hierarchy of 
terms and to a single scale governing the large 
time expansion of the energy density. 

In the absence of a scaling argument, one can 
still determine the small r behavior of the energy 
density e(r) assuming that at r = the initial 
condition is regular. This assumption does not 
need to be a priori relevant for realistic heavy-ion 
collisions, but it is in any case an interesting anal- 
ysis in its own. It implies in fact a finite limit of 
e(r) as r — >• 0, consistent with the s — behavior 
found in [22]. Moreover, it restricts the expansion 
of the energy density to be only in even powers 
of the proper-time 



The coefficients t% n are uniquely determined, 
through Einstein's equations, in terms of the co- 
efficients of the initial condition for the metric 



ao(z) = a(r = 0, z) 



(93) 



and similarly for the other components. 

The admissible initial conditions are given by 
the constraint equations contained in the full set 
of Einstein's equations. These imply ao(z) = 
bo(z). Defining 



«(* ) 
w{z 2 ) 



1 

4~ 



the remaining constraint reads 



W 



0. 



(94) 
(95) 

(96) 



where the prime denotes the z 2 -derivative. It is 
easy to see that there does not exist an anywhere 
bounded solution of the constraint equation. In- 
tegration of ([96]) gives 



= / (v'+w')dz z + / {v*+w*)dz 2 
Jo Jo 



(v 2 + w 2 )dz 2 



(97) 



since the first integral vanishes because of the 
boundary conditions. The only regular solution 
is vacuum AdSs with v = w = 0. Therefore, at 
any time, including r ~ 0, the metrics of interest 
are singular. 

The requirement imposed in |47] of the absence 
of curvature singularities (other than the one at 
z = oo) is thus a very powerful tool to select the 
allowed initial conditions. The coordinate singu- 
larity that is then left at all times in the bulk 
might signal the presence of a dynamical horizon, 
although this interesting statement requires fur- 
ther investigation. 

A bonus that follows from the simplicity of (|Mj) 
is that it can be solved analytically. One neat ex- 
ample of a solution to the initial value constraint, 
satisfying the non-singularity requirement, reads 



e(r) = e + c-it 2 + e 4 r 4 



(92) 



ao(z) = bo(z) = 21ncosaz 2 , 
Cq(z) = 2 In cosh az 2 , 



(98) 
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Figure 5. Pade approximation of the effective 
power s e ff(r) in the energy density for expan- 
sions up to order t 64 (dotted line), r 80 (dashed 
line), r 96 (solid line) and initial conditions ([98]). 
The two horizontal lines denote s = 1 (free 
streaming scenario) and s — 4/3 (perfect fluid 
case). 



for a generic constant a. The details of a larger 
number of allowed initial conditions can be found 
in [47]. Some of them are also mentioned in the 
following. 

Given the class of acceptable initial conditions, 
one then needs to solve Einstein's equations with 
these initial data to obtain the profile e(r). In the 
absence of a scaling variable, the solution has to 
be exact and the ideal treatment of the equations 
should be numerical. This analysis is in progress 

In a primary exam, an analytic analysis was 
pursued in [57J where the equations were solved 
in a power series in z and r for specific initial 
conditions. The limit of this approach is that the 
power series for the energy density has a finite 
radius of convergence necessitating the use of as 
a resummation scheme, such as Pade approxima- 
tion. 

The system evolves from an early time dynam- 
ics governed by the initial conditions to an hy- 
drodynamical regime at late times. One way to 
quantify this transition is through the exam of the 
effective exponent of the power-law dependence of 
the energy density 

s ef f(T) = -T— lne(r). (99) 




Figure 6. Pade resumed profiles for e(r) for the 
initial conditions v + w = tanh(z 2 ) — tan(z 2 ), 
v + w = tanh(z 2 + z 8 /6) — tan(z 2 ) and v + w = 
2/3z 6 (l + z 2 /2)/(z 2 - 1), respectively. 



In FigJS] this effective power is plotted as a func- 
tion of t for the initial conditions (|98|) . The plot 
is taken from [37]. The different curves corre- 
spond to different cuts in the Pade approxima- 
tion. All profiles start out at zero, clearly cross 
the line s = 1 denoting the free streaming sce- 
nario, and move upwards towards s = 4/3, which 
corresponds to the perfect fluid flow. However, to 
be sure that the profile indeed reaches s = 4/3, a 
numerical solutions is required. 

One can also perform the Pade approximation 
assuming the late time exponent s = 4/3. The 
profiles of e(r) obtained in this way are plotted 
in Fig E] for a set of initial conditions (plots taken 
from [47]). The energy densities differ in the ini- 
tial stages of the evolution, whereas in the late- 
time regime they seem to approach local equilib- 
rium. 

Recall that the positivity condition T^ff > 
implied -4c/t < e' < (see eq. {(HSJ))- In 
the third plot of FigE] a temporary violation of 
the positive energy condition is observed. Such 
a transient behavior may appear for a quantum- 
driven process (see e.g. the discussion in [49] , 
quoting field-theoretical results [50].). 

An alternative way to describe the transition to 
hydrodynamics is through the relative difference 
between the longitudinal and transverse pressures 
defined as 

Ap(r) = 1 - . (100) 

When the quantity Ap(r) approaches zero, it sig- 
nals isotropization indicating local equilibrium, 
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Figure 7. Relative difference in pressures 
Ap(r) = 1 — p\\/p± for the initial conditions 
v + w = tanh(z 2 ) — tan(2 2 ) and v + w = tanh(z 2 + 
z 8 /6) — tan(z 2 ), respectively. The profiles exhibit 
a rapid fall-off at r ~ 0(1), but the first profile 
does not reach yet complete isotropization. 



while a value of order one indicates a situation 
close to a free streaming scenario. FigJT] appears 
in [47] and shows the relative difference of pres- 
sures (|f OOj) for a couple of initial conditions. It is 
interesting to observe that both profiles exhibit 
a rapid fall-off on a scale r ~ 0(1)- In the 
first profile however, isotropization remains in- 
complete until r ~ 0(5). It is interesting to note 
that figures [5] and [7] arc consistent with the full 
anisotropy solution at very small proper-times 
but show a variety of behavior in the whole of 
the early proper-time region. 

These results point towards a fast isotropiza- 
tion, which could however remain incomplete for 
some time. The convergence to the hydrodynamic 
regime remains an open problem. In the coming 
more systematic numerical study |48j the stabil- 
ity of the Pade approximation has been shown 
reliable to a larger proper-time interval through 
non trivial reparamctrization of the metric. A 
more general proof of the transition to the hy- 
drodynamic regime, filling the gap between early- 
time and late-time behavior is still lacking. Some 
progress has been made, as we shall see now, 
starting from specific shock wave initial condi- 
tions mimicking the initial ultra-relativistic col- 
liding heavy-ions. 

4.3. Shock wave initial conditions 

An interesting way to model ultra- 
relativistically boosted large nuclei is through 
gauge theory shock waves. In this context, by 



shock wave we mean a plane shell of matter mov- 
ing at the speed of light. The gravity dual of a 
boundary shock wave has been first constructed 
in |15j and later applied to various analysis, 
such as in attempts to describe holographically 
deep inelastic scattering or heavy-ion reactions 
through shock waves collisions. 

In terms of the light-cone coordinates x = 
t ± y, an ultra relativistic source extended in 
the transverse direction (a model for an incident 
heavy ion with large transverse size) and moving 
along the direction x + has an energy-momentum 
tensor 

T__ = F(x~) ~ fj,5(x~). (101) 

The last term represents more specifically a shock 
wave, i.e. a shell with vanishing thickness and 
transverse energy density /i. It may be more 
concretely represented by a Gaussian centered at 
x- = 0. 

For arbitrary profile F(x~), the dual back- 
ground to this configuration reads 

z 2 ds 2 = -dx~ dx + +z i F(x~)dx~ 2 +dx 2 L +dz 2 (102) 

where F(x~) = (T ), and is an exact solution^"! 

to Einstein equations |15j . 

These solutions constitute the starting point 
towards solving the problem of shock waves colli- 
sions. Shock waves having a very small thickness 
along the collision axis and localized in the trans- 
verse directions qualitatively resemble the highly 
Lorentz-contractcd nuclei in a heavy-ion collision. 
FigH] appeared in [52] and depicts the space-time 
picture of ultra-relativistic heavy-ion collisions in 
the center of mass frame. In the AdS/CFT lan- 
guage, the collision of shock waves in Af = 4 SYM 
translates into colliding gravitational shock waves 
in the bulk. 

Several deserving attempts, both analytical 
and numerical, have been made to solve this prob- 
lem |49l52l53l54l55l56l57l58l59l60| . but the full 

understanding of the transition to the hydrody- 
namic regime is still an active subject of research. 
Let us discuss some of these angles of attack. 

10 Note that the solution (1 1 02 1) has been provided a gener- 
alization to profiles F = F(x~ , x± , z) in 1511 . 
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Figure 8. Space-time diagram of the collision of 
two ultra-relativistic nuclei in the center of mass 
frame. 



The collision of two singular shock waves mov- 
ing in the directions x ± on the boundary of a five- 
dimensional AdS space was studied in [49]. The 
pre-collision line element is given by the super- 
position of two metrics of the form (|102p . with 
delta-function Ansatz F(x ± ) ~ 6(x ± ). Solving 
Einstein's equations in a power series expansion 
in the proper-time allows to extract the expres- 
sion of the energy-momentum tensor right after 
the collision. In [35]: it was observed that the 
latter quantity is not boost-invariant in the sense 
of Bjorken. Boost-invariancc could be eventu- 
ally re-established at later times, at which the 
expansion breaks down. A crude estimate of the 
thermalization time in this simple model leads to 
Ttherm ~ (J-^ ■ The non-trivial dimcnsionlcss fac- 
tor in the above relation is 0(1). 

In [53] , it was observed that the metric of a sin- 
gle shock wave, moving along the light cone, can 
be computed exactly in first order perturbation 
theory. The only diagram contributing represents 
a single graviton exchange between the source nu- 
cleus at the AdS boundary and the location in 
the bulk where the metric is measured. The met- 
ric after the collision in the forward light cone 
can then be constructed perturbatively through 
an expansion in graviton exchanges. This com- 
putation was carried out at tree level in |52| and 
up to NNLO in [53], where the technique was ap- 



plied to proton- nucleus collisions. 

The analysis of the holographic stress-energy 
tensor leads to the strong prediction of full nu- 
clear stopping of the two nuclei shortly after the 
collision. This happens after a time of the order 
of the inverse typical transverse momentum scale, 
independently of the energy density. This behav- 
ior was interpreted in |52| as signaling black hole 
creation in the bulk. If the two nuclei stop com- 
pletely, the strong interactions exerted between 
them are likely to thermalize the system, even- 
tually leading to Landau hydrodynamics at late 
times [6]. This non boost-invariant solution is 
not expected to satisfy the phenomenological con- 
straints, even if the description of multiplicities 
seems to be valid. 

In [52|53j it was pointed out that weak inter- 
actions are expected to be determinant in the 
very early stages of the collision [61]. Since a 
realistic model should take them into account, 
and in absence of a dual description incorporating 
weak coupling aspects of QCD, it was proposed 
to mimic these effects through net zero-energy 
unphysical shock waves with profiles F(x ± ) ~ 
S'(x ). The resulting energy density of the 
strongly coupled medium starts out as a constant 
at early proper time, in agreement with the re- 
sults of [2"2] . 

The above analytic approaches have the limit 
of not being able to address the non-linearity of 
the thermalization process. A numerical treat- 
ment of the full problem appeared in [62], after 
these lectures were given|H] The main advantage 
of solving the initial value problem numerically 
is that it is more likely to properly describe the 
transition between the early and late time behav- 
iors. 

Ref. [62] considered the collision of two planar 
sheets in J\f = 4 SYM propagating towards each 
other at the speed of light, with finite energy den- 
sity, finite thickness and Gaussian profile 

F ( x ±) = ^^ e -(* ± /"V2)\ (103) 

V 2-KLU 2 

The equivalent gravitational problem is the col- 
lision of two planar shock waves which are reg- 

11 Other recent papers related to our lectures are quoted 
in Refs. |63l64l65l . 
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ular, non-singular, source-less solutions to Ein- 
stein's equations. The collision leads to horizon 
formation in the bulk and the numerical integra- 
tion was stopped in |62] at the apparent horizon 
location. 

The energy density of the collision products is 
found to be peaked around two receding maxima 
which move outwards at less than the speed of 
light. Setting the width in (fTUcif to oj = 0.75/ /x, 
the application of this model to the RHIC colli- 
sions leads to /i ~ 2.3 GeV. The total time re- 
quired for thcrmalization, from when the Gaus- 
sian shock waves start to overlap significantly un- 
til the onset of hydrodynamics, is estimated to be 
Autot ~ 4//x ~ 0.35 fm/c (v denotes time in in- 
going Eddington-Finkclstcin coordinates). This 
timescale is consistent with the results of [45] . 
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